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Abstract 

Motivated by concerns regarding possible two-body contributions to the recently-measured in- 
clusive T(nS)^ 7 + X direct photon spectrum, we report on a preliminary new study of exclusive 
radiative decays of the T resonances into two-body final states 7^7, with 1Z some resonant hadronic 
state decaying into four or more charged particles. Such two-body processes are not explicitly ad- 
dressed in the extant theoretical frameworks used to calculate the inclusive direct photon spectrum. 
Using data collected from the CLEO III detector at the Cornell Electron Storage Ring, we present 
upper limits for such T(1S), T(2S), and T(3S) two-body decays as a function of the recoil mass 
Mfc. Additionally, we place upper limits on the cross-section for 1Z production via radiative return 
for center-of-mass energies just below the T(1S), T(2S), and T(3S) energies. The results presented 
in this document are preliminary. 



* Submitted to the 33 rd International Conference on High Energy Physics, July 26 - August 2, 2006, Moscow 
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Introduction 

Very recently, we extracted a s from a measurement of the direct photon spectra in 
T(1S, 2S, 3S) — > gg-y To extrapolate beyond the experimentally accessible direct photon 
momentum region, those measurements relied on several theoretical parameterizations of 
the expected photon momentum spectrum in the T system 0, HI to obtain the total direct 
Y ~~ > 991 decay width relative to the dominant T — ► ggg width. The ratio of those widths 
can then be used to estimate the strong coupling constant at the energy scale of the T. These 
theoretical parameterizations generally neglect the contribution to the photon momentum 
spectrum due to two-body decays, e.g. T — > ggj — > TZj, with 1Z some resonant hadronic 
state. 

As noted in the a s extraction, contributions from such possible two-body decays may 
result in a slight underestimate of the extracted value of ot s . This systematic consideration 
in the gg^y analysis motivated a search for T — > TZj. We concern ourselves with high 
multiplicity (> 4 charged tracks) final states, as we employ the same hadronic event selection 
cuts in this analysis that we did in the ggj analysis We note that, although two-body 
branching fractions have been observed for, e.g., T(1S)— > 7/2(1270) at the level of 10 -4 , the 
fraction of ^(1270) decays into >4 charged tracks is only « 3% p. Our sensitivity to the 
/2(1270) is additionally reduced in this analysis by its large width (compared to our typical 
photon energy resolution). 

The analysis, in general terms, proceeds as follows. After selecting a high-quality sample 
of e + e~ annihilations into hadrons using the hadronic event selection cuts of the previous 
analysis pj, we plot the inclusive isolated photon spectrum in data taken at both on-T- 
resonance and off-T-resonance energies (the latter sample is used for systematic checks of 
the overall procedure). A two-body radiative decay of the T will produce a monochromatic 
photon in the lab frame; the momentum of the radiated photon is related to the mass of the 
recoil hadron 1Z via Equation [T] (x 7 = E 1 jE\ )eam ). 

M n = 2E heam ^l - (x^y (1) 

Assuming that the intrinsic width of the recoil hadron is much smaller than the exper- 
imental photon energy resolution, then the measured radiative photon energy should be a 
Gaussian centered at the momentum x?~. For a 1 GeV (4.5 GeV) recoil photon, this im- 
plies a recoil mass typically much narrower than 20 MeV (60 MeV). Not knowing a priori 
the mass of the hadron 7Z, we therefore perform a set of fits to the T(1S) photon spec- 
trum to a Gaussian signal, centered at the value x^, and with a resolution corresponding 
to the known CLEO-III electromagnetic calorimeter resolution, atop smooth polynomial 
backgrounds, over the 0.2 < x 7 < 1.0 photon energy range. 1 

We construct 95% confidence level upper limits from these fits by adding 1.645 * cr A (x^) 
to the A(x y ) distribution, where A{x 1 ) is the x 7 -dependent Gaussian fit area and a^ix^) is 
the fit error. We then recast this upper limit as a function of recoil mass M-r, correct for 



It should be noted that 'bumps' in the inclusive photon spectrum can be due not only to resonant two- 
body decays but also to continuum threshold effects like the crossing of the cc threshold. In the previous 
analysis, we also identified an excess of photons in data as x 7 — ► 1. Further examination of these events 
indicated that they were dominated by e + e~ — » 77r + 7r~7r + 7r~, although the possibility that the 4-pion 
state resulted from the decay of an intermediate resonance 1Z was not investigated. 
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the efficiency loss due to the fiducial acceptance of the detector (for the purposes of this 
correction, it is assumed that TZ is pseudo-scalar with a corresponding 1 + cos 2 9 angular 
distribution for the recoiling 7) and of the event and photon selection cuts that define 
our data sample. Note that the exact form of the efficiency correction due to the event and 
photon selection cuts varies with the decay final states considered for TZ. To be conservative, 
we derive our x 7 -dependent efficiency correction from the decay mode yielding the lowest 
reconstruction efficiency imaginable. 

This final efficiency-corrected limit is converted into an M^-dependent branching ratio 
upper limit Bi^ylZ) by dividing the resulting yield by the calculated total number of resonant 
T events. For the off-resonance running, the distributions are divided by the off-resonance 
luminosity for comparison's sake. An example of simulated signal superimposed on back- 
ground is given in Figure [T]for T(4S') — > 7 + TZ, TZ —* 7r + 7r~7r + 7r~. 




FIG. 1: T(45*) — > TZ^y, TZ — > 4ir, for various hypothetical TZ masses. The magnitude of B(T(4S) — > 
TZj, TZ — > 4-7r) has been grossly exaggerated for the sake of presentation. The red (dotted) curve 
is T(45) data, while the blue (solid) curve is T(45) data with signal Monte Carlo added on top. 
From the left, the first 'bump' corresponds to a TZ of mass ~ 7.6 GeV, the second from the left 
corresponds to a TZ of mass ~ 6.5 GeV, the third from the left corresponds to a TZ of mass ~ 4.7 
GeV and the fourth from the left corresponds a TZ of mass ~ 3.3 GeV. 



Event Selection 

Event selection criteria in this analysis are identical to those imposed in the previous 
analysis pj. The inclusive photon spectra are therefore identical to those taken from our 
previous analysis, as well. The background shape is approximately exponential in the region 
from 0.2 < x~ < 1.0. 
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Fitting the Inclusive Photon Spectrum 

To extract the possible magnitude of a two-body radiative signal, we step along the 
inclusive photon spectrum over the interval 0.2 < a; 7 < 1.0, fitting it to a Gaussian with 
width equal to the detector resolution at that value of photon energy, plus a background 
parametrized by a smooth Chebyshyev polynomial. We assume that the intrinsic width 
of the resonance 1Z is considerably smaller than the detector resolution. Our step size is 
determined by the energy resolution of the detector o"e; we use steps of width oe/2. For 
photons in the central "barrel" region of the Csl electromagnetic calorimeter, at energies 
greater than 2 GeV, the energy resolution is given by 

f(%) = ^ + 1.14-0.01£, (2) 

where E is the shower energy in GeV. At 100 MeV, the calorimetric performance is about 
20% poorer than indicated by this expression due to the material in front of the calorimeter 
itself. 

At each step, we use a ±10<7e fitting window; the background is expected to be relatively 
smooth over such a limited interval. We use a 3rd-(2nd-)order Chebyschev polynomial in the 
photon interval x 7 < (>)0.6. For each fit, the Gaussian fit area A(x^) and fit error ox(x 7 ) is 
recorded. Note that the fits are highly correlated point-to-point, and that the bin width is 
much finer than the detector resolution. Note also that we have not attempted to analytically 
correct for smearing-induced distortions that may result from the finite resolution of the 
detector; we simply assume that the smeared spectrum can be well- described by a smooth 
higher-order polynomial. 

Extracting Upper Limits 

To convert the A(x^) distribution obtained from fitting the inclusive photon spectrum 
into a one-sided upper 95% confidence interval upper limit, we add 1.645 *cr^(a; 7 ) point-wise 
to the v4(x 7 ) distribution, as a function of photon momentum. In this process, since we are 
interested in enhancements in the inclusive photon spectrum, all negative areas from the 
raw fits are set equal to 1.645 * a^x^) at these points. The resulting contour for the T(liS') 
fitting is shown in Figure 121 

We convert the limits as a function of photon energy x 7 into a function of a hypothetical 
resonance recoil mass M-r via Equation ^ For the purposes of this conversion the mean 
values for each running period of Eb eam are used for each data sample; we neglect the MeV- 
scale variation in beam energies for a particular run period. These values are given in 
Table I. 

The resulting M^-dependent contour, for the T(IS'), is shown in Figure EH 
Efficiency Correction 

We consider two efficiency corrections to the upper limit contour: one due to the fiducial 
acceptance of the detector, and the other due to our event and shower selection cuts. 

For photons in the "barrel" of the detector (| cos#| < 0.707, with 9 the polar angle of 
the photon momentum relative to the beam axis), we assume that TZ is a pseudo-scalar, 
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FIG. 2: A(xj) + 1.645 * <ta(x 7 ) versus x 7 for fits to the Y(1S) inclusive photon spectrum, where 
negative points have been mapped to 1.645 * 07i(x 7 ). 



Event Type 


Average Beam Energy (E beam ) 


Resonance Y(1S) 


4.73 


GeV 


Resonance T(2S) 


5.01 


GeV 


Resonance T(3S) 


5.18 


GeV 


Continuum T(IS') 


4.72 


GeV 


Continuum T(2S') 


5.00 


GeV 


Continuum Y(3»S) 


5.16 


GeV 


Continuum T(45) 


5.27 


GeV 



TABLE I: The mean values of Eb eam for each data sample used in this analysis, used to map our 
upper limit contour from a function of x 7 to a function of M-ji. 

which corresponds to a 1 + cos 2 9 distribution on the photons in the two-body decays we are 
considering. This choice is arbitrary and amounts to an 0.6 uniform acceptance efficiency 
correction to our limit. In addition to this angular acceptance correction, we assess an 
efficiency correction due to the CLEO III detection efficiency. Not knowing a priori what 
the decay mode of our hypothetical resonance 1Z will be, we have generated 5000-event Monte 
Carlo samples spanning a wide range of final state charged (and neutral) multiplicities and 
masses M n . In the interests of producing a conservative upper limit, we take our efficiency 
correction from the mode which yields the lowest average efficiency in the main M-ji interval 
of interest. A list of 1Z decay modes considered in this analysis and their average efficiencies 
(averaged over the photon momentum spectrum from 1.0 GeV < E 1 < 4.5 GeV) is given in 
Table II. We thereby derive a collection of efficiency measurements distributed over photon 
momentum E 7 , each of which corresponds to a different hypothesis for the mass M n . To 
obtain efficiencies between these points (where we have not explicitly filled in the efficiency 
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FIG. 3: A(Mn) + 1.645 *(ta(Mii) versus M-ji for fits to the Y(IS') inclusive photon spectrum. This 
curve was obtained from mapping the points in Figure [21 into recoil mass space via Equation ^ 

with a M-ji hypothesis), we perform a linear extrapolation between the two neighboring 
points. In this manner, we point-wise efficiency correct our upper limit as a function of 
x 7 (or Ey) before mapping the upper limit into Mn. We find that the lowest efficiency for 
various hypothetical decay modes of 1Z considered was obtained from 1Z — ► 4Kn° (Figure 
We therefore use this efficiency function to determine our upper limit contours. 

Results 

To convert the efficiency corrected upper limit contour into an upper limit on the two- 
body radiative branching ratio B(jlZ), we simply divide the efficiency-corrected upper limit 
contour by the total calculated number of T(1S), T(2S) and T(3S) events. These numbers 
were calculated in our previous analysis ^ , and for convenience are reproduced in Table III. 
For comparison's sake, the efficiency-corrected upper limit contours obtained from fitting 
below T resonance data have been divided by the total luminosity of off-resonance data 
taking and therefore correspond to an upper limit on a cross-section. The resulting on- 
resonance upper limit Bi^ylZ) is shown in Figure and the resulting off-resonance cross- 
section "limit" is shown in Figure |3 To set the scale of the continuum cross-section limits, 
the raw ISR cross-section for e + e~ — > ip + 7 is expected to be ~5 pb in the 10 GeV center- 
of-mass regime, implying an expected observed cross-section into >4 charged tracks ~ 10 -4 
nb, taking into account the strong forward peaking expected for ISR processes. 

Given the fact that we have not performed a continuum subtraction on the on-resonance 
inclusive photon spectrum from T decays, it is interesting to compare the structure observed 
in Figure with structure observed when we apply the fitting procedure to continuum data. 
Figures |H1 and M show the resonances limits of Figure El separately, with the respective 
continuum limits of Figure El overlaid for comparison. We observe some, albeit not complete, 
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Event Type 


Average Efficiency (e) 


n -» K+K--K+-K- 


0.53 


± 


0.03 




0.53 


± 


0.02 


71 -»■ K + ir-7r+vr-7r 7r 


0.54 


± 


0.02 


72. -> K+K-p+p- 


0.56 


± 


0.02 


1Z -» K + K-p + p--K° 


0.50 


± 


0.05 


TZ -► K + K-p+p-ir°ir 


0.57 


± 


0.02 


IZ — > p + p _ 7r + 7r _ 


0.62 


± 


0.03 


72. — > p + p~TT + TT^TT° 


0.54 


± 


0.05 


IZ — > p + p _ 7r + 7r _ 7r 7r 


0.63 


± 


0.02 


72. -» K+K-K+K- 


0.50 


± 


0.02 


IZ -> K + K-K + K--k°tt° 


0.49 


± 


0.02 


IZ — ► p + p~p + p~ 


0.67 


± 


0.02 


TZ — » p + p~p + p~7T° 


0.65 


± 


0.02 


72 — > p + p~p + p~ir°ir 


0.63 


± 


0.02 


IZ — > 7r + 7r _ 7r + 7r _ 


0.59 


± 


0.02 


72. — > 7T + 7r _ 7r + 7r _ 7r 


0.65 


± 


0.02 


72. — > 7T + 7r _ 7r + 7r _ 7r 7r 


0.59 


± 


0.01 


72. — > 7r + 7r _ 7r + 7r~47r 


0.57 


± 


0.02 


72 — > 7r + 7r~7r + 7r"~67r 


0.60 


± 


0.02 


72. — > 7T + 7r _ 7r + 7r _ 87r° 


0.60 


± 


0.02 


72. -» K + K-K + K~K + K- 


0.68 


± 


0.04 


1Z — > p + p~p + p~p + p~ 


0.53 


± 


0.04 


IZ — > 7r + 7r _ 7r + 7r _ 7r + 7r _ 


0.74 


± 


0.03 


72- -» K+K-K+K-tt 


0.48 


± 0.02 



TABLE II: Average efficiencies for the reconstruction of various events that could be detected in 
this analysis, obtained by fitting the photon momentum-dependent reconstruction efficiencies to a 
straight line in the interval 1.0 GeV < E~ < 4.5 GeV. The presented errors on the efficiencies are 
statistical only. The lowest-efficiency final state (K + K~K + K~ir°) is used for setting upper limits. 

correlation between the two. 



T Resonance 


AWT(nS)) (xlO 6 ) 


T(1S) 
T(2S) 
T(3S) 


21.0 ± 0.06 
8.4 ± 0.04 
5.2 ± 0.06 



TABLE III: The total number of calculated T(1S), T(2S) and T(3S) events in our data samples 
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FIG. 4: The efficiency for detecting an T -> 7 + ft, ft -> 4K7T event as a function of observed 
photon energy E~. Each point in this efficiency is obtained from a different ft mass hypothesis. 
This photon momentum-dependent efficiency correction distribution is used to point-wise correct 
our upper limit, where the efficiency between points in this distribution is estimated by linear 
interpolation. 




FIG. 5: The M^-dependent B{^1Z) upper limit contours obtained for T — > 7 + ft, ft — >> 4 charged 
tracks for the T(LS), T{2S) and T(35). All limits are of order B(-yK) fa 1(T 4 . 
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FIG. 6: Our M^-dependent cross-section upper limit contours obtained for e + e~ — > j+TZ, 1Z — >> 4 
charged tracks for the continua below T(15), T(25), T(3S) and T(4S*) (nb). This plot is obtained 
by dividing the result of our fitting procedure on the continuum by the off-resonance luminosity. 
The angular correction here is (as before) 1 + cos 2 6L. 
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FIG. 7: Comparison of the M^-dependent Gaussian fit area A(Mn) for the T(15) versus the 
continuum below T(1S). Note the correlation between the resonance and the below-resonance 
structure. Note also that the below-resonance upper limit curve has been scaled to match the 
'floor' of the on-resonant upper limit curve. 
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FIG. 8: Comparison of the M^-dependent Gaussian fit area A(Mn) for the T(2S) versus the 
continuum below T(25). Note that the below-resonance upper limit curve has been scaled to 
match the 'floor' of the on-resonant upper limit curve. 
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1000 
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4.2 6.2 

M H (GeV) 



FIG. 9: Comparison of the M^-dependent Gaussian fit area A{Mn) for the T(35) versus the 
continuum below T(35). Note that the below-resonance upper limit curve has been scaled to 
match the 'floor' of the on-resonant upper limit curve. 
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Cross-Check 

In order to ensure that we are able to identify a signal at a given sensitivity level, we 
embedded pure Monte Carlo signal in data, and performed our fitting procedure on the 
resulting distribution in order to ensure that we recovered the correct signal magnitude in 
our branching ratio upper limit. To do this, T(AS) — > 7 + 7Z, TZ — > 7r + 7r~7r + 7r _ events were 
embedded in the T(4S I ) inclusive photon spectrum with branching ratios of the order of 10~ 5 , 
1(T 4 , l(r 3 , and 10" 2 , under 10 different M n hypotheses: M n = 0.6 GeV, 1.5 GeV, 2.5 GeV, 
3.5 GeV, 4.5 GeV, 6.5 GeV, 7.5 GeV, 8.5 GeV, 9.5 GeV and 10.5 GeV. The resulting upper 
limit contours derived from applying our procedure to these spectra are show in Figure [TDJ 
We reconstruct all signals above our upper limit floor (around 10 -4 ) that are within our 
accessible recoil mass range. 




FIG. 10: Upper limit contours derived from applying our procedure to fabricated Monte Carlo 
signal spectra. We reconstruct all input signals above our upper limit floor (~ 10 -4 ) that are 
within our accessible recoil mass range. 



Systematic errors 

We identify and account for systematic errors as follows: 

1. We account for possible systematics in our event and shower reconstruction efficiency 
by using the lowest-efficiency final state considered. 

2. We assess a uniform 1% degradation of the limit due to the luminosity uncertainty as 
calculated in the previous analysis [l|. 

3. We degrade the limit uniformly by the uncertainty in the calculated number of total 
T(1S), T{2S) and T(35) events. 
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Summary 



As shown in Figured! our sensitivity is of order 10 -4 across all possible values of Mr, well 
above the threshold for any known branching ratio for T — > 7+pseudoscalar, pseudoscalar^ 
h + h~ h + h~ +neutrdls. We measure upper limits of &(T(1S) — > 7 + 71,71 — >> 4 charged 
tracks ) < 1.05 x 1(T 3 , B{T(2S) -> 7 + 71,71 ^> 4 charged tracks ) < 1.65 x 10~ 3 and 
B(T(3S) — > J+7Z, 7Z ^> 4 charged tracks ) < 2.47 x 1CT 3 for all possible masses Mr, under 
the assumption that 7Z is a pseudoscalar. Constraining 1.5 GeV < Mr < 5.0 GeV we set a 
much more stringent limit of B(T(1S) — > 7 + 7Z, 7Z — >> 4 charged tracks ) < 1.82 x 10~ 4 , 
B(T(2S) -> j + 71,71 -►> 4 charged tracks) < 1.69 x 10~ 4 and B(T(3S) -> 7 + ft,ft^>4 
charged tracks ) < 3.00 x 10~ 4 . Additionally, we report these upper limits as a function of 
the mass recoiling against the photon, as shown in Figure 03 We also report cross-section 
limits for resonant processes on the continuum (Figure 

We note that we appropriately ignored the distortion of the inclusive photon spectrum 
in our previous extraction of a s due to two-body decays as we limit the branching ratio of 
these events to be < 10~ 4 . Further work on exclusive multiparticle final states (e.g., 7471", 
jiK, ^K°K° and •yK°K7i) would help elucidate the nature of such radiative decays. 
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